Comparison with Kapsenberg et al. (2017)
Description
Time series trends analysis by residuals (substracting climatological monthly means)
Kapsenberg et al. (2017)
2007-2015 : Comparison with Kapsenberg results
SURFACE
Table of monthly means (2007-2015) - surface data
➝ Comparaison supplement, Tab S1 (Kapsenberg) same period : results ok ✓
monthly_means_0m_0715 <- ungroup(data_0715) %>%
group_by(monthd) %>%
summarise(
Salinity_month = mean(salinity, na.rm = TRUE),
Temperature_month = mean(temperature, na.rm = TRUE),
dic_month = mean(dic, na.rm = TRUE),
ta_month = mean(ta, na.rm = TRUE),
pH_calc_month = mean(pH_calc, na.rm = TRUE),
pH18_month = mean(pH18, na.rm = TRUE),
pCO2_month = mean(pCO2, na.rm = TRUE),
OmegaCalcite_month = mean(OmegaCalcite, na.rm = TRUE),
OmegaAragonite_month = mean(OmegaAragonite, na.rm = TRUE),
Nta_month = mean(Nta, na.rm = TRUE),
Ndic_month = mean(dic, na.rm = TRUE))
monthly_means_0m_0715## # A tibble: 12 × 12
## monthd Salinity_month Temperature_month dic_month ta_month pH_calc_month
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 1 37.9 14.4 2251. 2551. 8.13
## 2 2 37.9 13.6 2260. 2556. 8.14
## 3 3 37.9 13.7 2260. 2555. 8.13
## 4 4 37.8 15.0 2254. 2553. 8.12
## 5 5 37.6 18.1 2241. 2548. 8.09
## 6 6 37.7 21.3 2233. 2544. 8.05
## 7 7 37.9 24.1 2228. 2549. 8.02
## 8 8 38.1 24.8 2227. 2556. 8.02
## 9 9 38.2 23.5 2230. 2562. 8.04
## 10 10 38.2 20.8 2227. 2561. 8.08
## 11 11 38.0 18.2 2232. 2555. 8.11
## 12 12 38.0 16.1 2242. 2556. 8.12
## # ℹ 6 more variables: pH18_month <dbl>, pCO2_month <dbl>,
## # OmegaCalcite_month <dbl>, OmegaAragonite_month <dbl>, Nta_month <dbl>,
## # Ndic_month <dbl>Table of time series anomalies regression analyses (2007-2015) - surface data
➝ Comparison Table 2 (Kapsenberg) same period : results ok ✓
ano_0m_0715 <- left_join(ungroup(data_0715), monthly_means_0m_0715, by = "monthd") %>%
mutate(
Salinity_ano = salinity - Salinity_month,
Temperature_ano = temperature - Temperature_month,
dic_ano = dic - dic_month,
ta_ano = ta - ta_month,
pH_calc_ano = pH_calc - pH_calc_month,
pH18_ano = pH18 - pH18_month,
pCO2_ano = pCO2 - pCO2_month,
OmegaCalcite_ano = OmegaCalcite - OmegaCalcite_month,
OmegaAragonite_ano = OmegaAragonite - OmegaAragonite_month,
Nta_ano = Nta - Nta_month,
Ndic_ano = dic - Ndic_month)
# Regressions and tables of key parameters
var_list <-
c("salinity",
"temperature",
"dic",
"ta",
"pH_calc",
"pH18",
"pCO2",
"OmegaCalcite",
"OmegaAragonite",
"Nta",
"Ndic")
lm.test <- vector("list", length(var_list))
for(i in seq_along(var_list)){
lm.test[[i]] <- lm(reformulate(var_list[i], "sampling_date"), data = ano_0m_0715)
}
lms <- lapply(var_list, function(x) {
summary(lm(substitute(i ~ decimal_date(sampling_date), list(i = as.name(x))),
data = ano_0m_0715))
})
reg <- NULL
for (i in 1:length(var_list)) {
#one loops through all anomalies
# calculate probability of fstatistic because it cannot be extracted from lms above
# see http://stats.stackexchange.com/questions/92824/extracting-the-model-p-value-for-a-multiple-regression-in-r
# slope returned by lm is per second
prob <-
pf(lms[[i]]$fstatistic[1],
lms[[i]]$fstatistic[2],
lms[[i]]$fstatistic[3],
lower.tail = FALSE)
reg <-
rbind(reg, as.numeric(
c(
lms[[i]]$coefficients[2, 1],
lms[[i]]$coefficients[2, 2],
lms[[i]]$coefficients[2, 4],
lms[[i]]$coefficients[1, 1],
lms[[i]]$coefficients[1, 2],
lms[[i]]$coefficients[1, 4],
lms[[i]]$fstatistic[1],
lms[[i]]$fstatistic[3],
lms[[i]]$r.squared,
prob
)
))
}
colnames(reg) <-
c("Slope",
"SE Slope",
"P Slope",
"Intercept",
"SE int.",
"P int.",
"F",
"df",
"R2",
"P value")
row.names(reg) <- var_list
reg_ano_surf <- reg
slope_temp <- reg_ano_surf[2, 1]
# Regression and table anomalies
var_list <-
c(
"Salinity_ano",
"Temperature_ano",
"dic_ano",
"ta_ano",
"pH_calc_ano",
"pH18_ano",
"pCO2_ano",
"OmegaCalcite_ano",
"OmegaAragonite_ano",
"Nta_ano",
"Ndic_ano"
)
lms <- lapply(var_list, function(x) {
summary(lm(substitute(i ~ decimal_date(sampling_date), list(i = as.name(x))),
data = ano_0m_0715))
})
reg <- NULL
for (i in 1:length(var_list)) {
#one loops through all anomalies
# calculate probability of fstatistic because it cannot be extracted from lms above
# see http://stats.stackexchange.com/questions/92824/extracting-the-model-p-value-for-a-multiple-regression-in-r
# slope returned by lm is per second
prob <-
pf(lms[[i]]$fstatistic[1],
lms[[i]]$fstatistic[2],
lms[[i]]$fstatistic[3],
lower.tail = FALSE)
reg <-
rbind(reg, as.numeric(
c(
lms[[i]]$coefficients[2, 1],
lms[[i]]$coefficients[2, 2],
lms[[i]]$coefficients[2, 4],
lms[[i]]$coefficients[1, 1],
lms[[i]]$coefficients[1, 2],
lms[[i]]$coefficients[1, 4],
lms[[i]]$fstatistic[1],
lms[[i]]$fstatistic[3],
lms[[i]]$r.squared,
prob
)
))
}
colnames(reg) <-
c("Slope",
"SE Slope",
"P Slope",
"Intercept",
"SE int.",
"P int.",
"F",
"df",
"R2",
"P value")
row.names(reg) <- var_list
reg_ano_surf <- reg
slope_temp <- reg_ano_surf[2, 1]
slope_pH <- reg_ano_surf[7, 1]
reg_ano_surf## Slope SE Slope P Slope Intercept
## Salinity_ano -0.001689114 0.0043739500 6.995639e-01 3.397891
## Temperature_ano 0.071869400 0.0216700091 9.924272e-04 -144.576487
## dic_ano 3.159510272 0.1933457350 1.541892e-46 -6355.699197
## ta_ano 2.080353058 0.1867599667 2.197247e-25 -4184.893358
## pH_calc_ano -0.002908521 0.0003337822 7.606349e-17 5.850844
## pH18_ano -0.001884556 0.0002215758 3.540344e-16 3.791014
## pCO2_ano 3.627942537 0.3896601695 8.181341e-19 -7298.048505
## OmegaCalcite_ano -0.012171909 0.0022940518 1.850430e-07 24.485279
## OmegaAragonite_ano -0.007271327 0.0015854405 6.022437e-06 14.627159
## Nta_ano 2.213010149 0.2738614887 7.131981e-15 -4451.749878
## Ndic_ano 3.159510272 0.1933457350 1.541892e-46 -6355.699197
## SE int. P int. F df R2
## Salinity_ano 8.7988267 6.995641e-01 0.1491318 416 0.0003583613
## Temperature_ano 43.5926333 9.924367e-04 10.9994216 411 0.0260650158
## dic_ano 388.9363365 1.542099e-46 267.0363299 408 0.3955880863
## ta_ano 375.6916128 2.197423e-25 124.0813516 414 0.2305996133
## pH_calc_ano 0.6714441 7.606755e-17 75.9308053 405 0.1578830146
## pH18_ano 0.4457270 3.540525e-16 72.3391874 405 0.1515467185
## pCO2_ano 783.8495602 8.181828e-19 86.6859612 405 0.1763035108
## OmegaCalcite_ano 4.6147686 1.850472e-07 28.1521109 405 0.0649935904
## OmegaAragonite_ano 3.1893094 6.022540e-06 21.0342677 405 0.0493722437
## Nta_ano 550.9074895 7.132321e-15 65.2987521 414 0.1362381017
## Ndic_ano 388.9363365 1.542099e-46 267.0363299 408 0.3955880863
## P value
## Salinity_ano 6.995639e-01
## Temperature_ano 9.924272e-04
## dic_ano 1.541892e-46
## ta_ano 2.197247e-25
## pH_calc_ano 7.606349e-17
## pH18_ano 3.540344e-16
## pCO2_ano 8.181341e-19
## OmegaCalcite_ano 1.850430e-07
## OmegaAragonite_ano 6.022437e-06
## Nta_ano 7.131981e-15
## Ndic_ano 1.541892e-46Plot of anomalies : temperature (2007-2015) - surface
Plot of anomalies : salinity (2007-2015) - surface
Plot of anomalies : pH (2007-2015) - surface
Plot of anomalies : alkalinity (2007-2015) - surface
Plot of anomalies : pCO2 (2007-2015) - surface
DEPTH 50 m
Table of monthly means (2007-2015) - 50m data
➝ Comparaison supplement, Tab S1 (Kapsenberg) same period : results ok ✓
## # A tibble: 12 × 12
## monthd Salinity_month Temperature_month dic_month ta_month pH_calc_month
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 1 38.0 14.5 2253. 2555. 8.13
## 2 2 38.0 13.7 2260. 2556. 8.13
## 3 3 38.0 13.6 2261. 2555. 8.13
## 4 4 38.0 14.1 2258. 2553. 8.12
## 5 5 38.0 15.0 2253. 2551. 8.12
## 6 6 38.0 15.1 2249. 2549. 8.12
## 7 7 38.0 15.4 2246. 2546. 8.11
## 8 8 38.0 15.3 2241. 2544. 8.12
## 9 9 38.0 16.2 2237. 2548. 8.12
## 10 10 38.1 18.5 2233. 2555. 8.10
## 11 11 38.1 18.0 2234. 2557. 8.11
## 12 12 38.0 16.2 2242. 2557. 8.12
## # ℹ 6 more variables: pH18_month <dbl>, pCO2_month <dbl>,
## # OmegaCalcite_month <dbl>, OmegaAragonite_month <dbl>, Nta_month <dbl>,
## # Ndic_month <dbl>Table of time series anomalies regression analyses (2007-2015) - 50m data
➝ Comparison Table 2 (Kapsenberg) same period: results ok ✓
## Slope SE Slope P Slope Intercept
## Salinity_ano 0.006280090 0.0020205419 2.013976e-03 -12.633358
## Temperature_ano 0.087907062 0.0187730250 3.867940e-06 -176.840381
## dic_ano 2.237743420 0.2042325862 1.192951e-24 -4501.572889
## ta_ano 1.612015807 0.1533532351 4.883020e-23 -3242.803820
## pH_calc_ano -0.002647707 0.0002445368 3.749752e-24 5.326314
## pH18_ano -0.001305946 0.0002666917 1.412206e-06 2.627133
## pCO2_ano 2.813812116 0.2453557568 1.528315e-26 -5660.463142
## OmegaCalcite_ano -0.007213340 0.0026443975 6.654962e-03 14.510865
## OmegaAragonite_ano -0.003929317 0.0018128827 3.078619e-02 7.904491
## Nta_ano 1.187830184 0.1244390877 1.254187e-19 -2389.492859
## Ndic_ano 2.237743420 0.2042325862 1.192951e-24 -4501.572889
## SE int. P int. F df R2
## Salinity_ano 4.0646310 2.013993e-03 9.660421 410 0.02301962
## Temperature_ano 37.7652440 3.868010e-06 21.926990 406 0.05124003
## dic_ano 410.8463170 1.193045e-24 120.052305 406 0.22821363
## ta_ano 308.4925535 4.883378e-23 110.497565 409 0.21270083
## pH_calc_ano 0.4919279 3.750038e-24 117.233404 403 0.22534771
## pH18_ano 0.5364963 1.412234e-06 23.979014 403 0.05615970
## pCO2_ano 493.5753768 1.528442e-26 131.521782 403 0.24605505
## OmegaCalcite_ano 5.3196612 6.655007e-03 7.440796 403 0.01812879
## OmegaAragonite_ano 3.6469260 3.078633e-02 4.697802 403 0.01152275
## Nta_ano 250.3274998 1.254266e-19 91.116090 409 0.18218988
## Ndic_ano 410.8463170 1.193045e-24 120.052305 406 0.22821363
## P value
## Salinity_ano 2.013976e-03
## Temperature_ano 3.867940e-06
## dic_ano 1.192951e-24
## ta_ano 4.883020e-23
## pH_calc_ano 3.749752e-24
## pH18_ano 1.412206e-06
## pCO2_ano 1.528315e-26
## OmegaCalcite_ano 6.654962e-03
## OmegaAragonite_ano 3.078619e-02
## Nta_ano 1.254187e-19
## Ndic_ano 1.192951e-24Plot of anomalies : temperature (2007-2015) - 50m
Plot of anomalies : salinity (2007-2015) - 50m
Plot of anomalies : pH (2007-2015) - 50m
Plot of anomalies : alcalinity (2007-2015) - 50m
Plot of anomalies : pCO2 (2007-2015) - 50m
2007-2022
→ SOMLIT clean + interpolations ok ☑
SURFACE
Table of monthly means (2007-2022) - surface data
Table of time series anomaly regression analyses (2007-2022) - surface data
## Slope SE Slope P Slope Intercept
## Temperature_ano 5.986065e-02 0.0090266801 5.975146e-11 -120.617872
## Salinity_ano 5.640707e-05 0.0020994838 9.785722e-01 -0.113659
## pH_calc_ano -3.964373e-03 0.0001547736 9.561431e-107 7.987629
## pH18_ano -3.113001e-03 0.0001180322 1.869102e-111 6.272240
## ta_ano 5.533254e-01 0.0961513994 1.226052e-08 -1114.869198
## dic_ano 2.439697e+00 0.1059150698 6.601288e-91 -4915.630099
## pCO2_ano 5.081188e+00 0.2104165592 6.711138e-95 -10238.096229
## SE int. P int. F df R2
## Temperature_ano 18.1886073 5.975814e-11 4.397711e+01 832 5.020349e-02
## Salinity_ano 4.2304243 9.785722e-01 7.218418e-04 832 8.675975e-07
## pH_calc_ano 0.3118469 9.570277e-107 6.560775e+02 819 4.447749e-01
## pH18_ano 0.2378181 1.870888e-111 6.955968e+02 819 4.592620e-01
## ta_ano 193.7314049 1.226154e-08 3.311692e+01 819 3.886429e-02
## dic_ano 213.4038134 6.606689e-91 5.305866e+02 819 3.931475e-01
## pCO2_ano 423.9880107 6.810126e-95 5.831370e+02 722 4.468014e-01
## P value
## Temperature_ano 5.975146e-11
## Salinity_ano 9.785722e-01
## pH_calc_ano 9.561431e-107
## pH18_ano 1.869102e-111
## ta_ano 1.226052e-08
## dic_ano 6.601288e-91
## pCO2_ano 6.711138e-95Plot of anomalies : Temperature (2007-2022) - surface
Plot of anomalies : salinity (2007-2022) - surface
Plot of anomalies : pH_calc (2007-2022) - surface
Plot of anomalies : pH(18) (2007-2022) - surface
Plot of anomalies : alkalinity (2007-2022) - surface
Plot of anomalies : DIC (2007-2022) - surface
Plot of anomalies : pCO2 (2007-2022) - surface
DEPTH 50 m
Table of monthly means (2007-2022) - 50m data
(rajouter les SD)
## # A tibble: 12 × 15
## Month Salinity_month SD_sal Temperature_month SD_temp dic_month SD_dic
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 1 38.1 0.123 14.6 0.863 2259. 18.4
## 2 2 38.1 0.121 13.8 0.388 2271. 17.1
## 3 3 38.1 0.134 13.7 0.387 2271. 15.7
## 4 4 38.0 0.113 14.2 0.521 2268. 17.6
## 5 5 38.0 0.120 15.1 0.806 2263. 19.7
## 6 6 38.0 0.0975 15.2 0.672 2258. 18.1
## 7 7 38.0 0.0826 15.4 0.920 2256. 19.0
## 8 8 38.0 0.0743 15.4 0.590 2253. 16.9
## 9 9 38.0 0.103 16.2 1.71 2246. 18.4
## 10 10 38.1 0.170 18.5 1.76 2236. 14.9
## 11 11 38.1 0.124 18.1 1.13 2237. 14.9
## 12 12 38.1 0.140 16.1 0.968 2249. 14.2
## # ℹ 8 more variables: ta_month <dbl>, SD_ta <dbl>, pH_calc_month <dbl>,
## # SD_PH <dbl>, pH18_month <dbl>, SD_PH18 <dbl>, pCO2_month <dbl>,
## # SD_pCO2 <dbl>Table of time series anomaly regression analyses (2007-2022) - 50m data
## Slope SE Slope P Slope Intercept
## Temperature_ano 0.026456564 0.0073678259 3.489534e-04 -53.309392
## Salinity_ano 0.009666510 0.0008266589 2.362696e-29 -19.477804
## pH_calc_ano -0.003343820 0.0001278701 4.579139e-110 6.737306
## pH18_ano -0.002980286 0.0001380500 3.547310e-82 6.004839
## ta_ano 0.583934397 0.0803203798 8.396021e-13 -1176.541935
## dic_ano 2.289966337 0.1044826549 3.779468e-84 -4613.945403
## pCO2_ano 3.454482642 0.1256618844 2.184751e-118 -6960.274502
## SE int. P int. F df R2
## Temperature_ano 14.8460452 3.489658e-04 12.89401 832 0.01526110
## Salinity_ano 1.6657037 2.363428e-29 136.73722 832 0.14114996
## pH_calc_ano 0.2576402 4.583474e-110 683.82918 819 0.45502788
## pH18_ano 0.2781511 3.549988e-82 466.06148 819 0.36267641
## ta_ano 161.8341636 8.397098e-13 52.85385 819 0.06062237
## dic_ano 210.5177179 3.782376e-84 480.36328 819 0.36969128
## pCO2_ano 253.1908589 2.186932e-118 755.71648 819 0.47990638
## P value
## Temperature_ano 3.489534e-04
## Salinity_ano 2.362696e-29
## pH_calc_ano 4.579139e-110
## pH18_ano 3.547310e-82
## ta_ano 8.396021e-13
## dic_ano 3.779468e-84
## pCO2_ano 2.184751e-118Plot of anomalies : Temperature (2007-2022) - 50m
Plot of anomalies : salinity (2007-2022) - 50m
Plot of anomalies : pH_calc (2007-2022) - 50m
Plot of anomalies : pH(18) (2007-2022) - 50m
Plot of anomalies : alkalinity (2007-2022) - 50m
Plot of anomalies : DIC (2007-2022) - 50m
Plot of anomalies : pCO2 (2007-2022) - 50m
1995-2022 : temperature & salinity
05-92 to 12-22 but not enough data before Jan. 1995
SURFACE
Table of monthly means for salinity and temperature (1995-2022) - 0m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.0 0.288 14.3 0.699
## 2 02 38.0 0.192 13.6 0.419
## 3 03 38.0 0.201 13.6 0.498
## 4 04 37.9 0.238 14.8 0.934
## 5 05 37.7 0.255 17.8 1.54
## 6 06 37.7 0.277 21.6 1.88
## 7 07 37.9 0.302 24.4 1.94
## 8 08 38.1 0.331 25.3 1.33
## 9 09 38.2 0.176 23.4 1.39
## 10 10 38.1 0.214 20.8 1.16
## 11 11 38.0 0.370 18.1 1.07
## 12 12 38.0 0.329 16.0 0.951Table of salinity and temperature time series anomaly regression analyses (1995-2022) - 0m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano -0.001033649 0.0009029316 2.525074e-01 2.076398 1.813827
## Temperature_ano 0.038182803 0.0040684385 2.570245e-20 -76.701761 8.172761
## P int. F df R2 P value
## Salinity_ano 2.525113e-01 1.310498 1350 0.0009697979 2.525074e-01
## Temperature_ano 2.571986e-20 88.080569 1350 0.0612487026 2.570245e-20Plot of anomalies : Temperature (1995-2022) - 0m
Plot of anomalies : salinity (1995-2022) - 0m
DEPTH 50
Table of monthly means for salinity and temperature (1995-2022) - 50m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.1 0.130 14.3 0.668
## 2 02 38.1 0.121 13.6 0.417
## 3 03 38.1 0.125 13.5 0.421
## 4 04 38.0 0.107 14.0 0.539
## 5 05 38.0 0.117 14.8 0.765
## 6 06 38.0 0.104 15.1 0.635
## 7 07 38.0 0.0982 15.2 0.633
## 8 08 38.0 0.101 15.4 0.689
## 9 09 38.0 0.110 16.1 1.54
## 10 10 38.1 0.165 18.6 1.81
## 11 11 38.1 0.129 17.9 1.14
## 12 12 38.1 0.132 16.0 0.953Table of salinity and temperature time series anomaly regression analyses (1995-2022) - 50m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano 0.003886329 0.0003942302 3.490107e-22 -7.806559 0.7919055
## Temperature_ano 0.019299605 0.0031399001 1.040676e-09 -38.767620 6.3072498
## P int. F df R2 P value
## Salinity_ano 3.49267e-22 97.18058 1347 0.06729115 3.490107e-22
## Temperature_ano 1.04099e-09 37.78028 1349 0.02724316 1.040676e-09Plot of anomalies : Temperature (1995-2022) - 50m
Plot of anomalies : salinity (1995-2022) - 50m
2000-2022 : Analysis every 10 years
SURFACE
Jan.2000 - Dec.2009
Table of monthly means (2000-2009) - 0m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.1 0.123 14.1 0.643
## 2 02 38.1 0.131 13.4 0.460
## 3 03 38.1 0.107 13.6 0.525
## 4 04 37.9 0.156 14.5 0.939
## 5 05 37.7 0.290 18.0 1.50
## 6 06 37.8 0.284 21.7 1.87
## 7 07 38.0 0.217 23.9 1.83
## 8 08 38.2 0.193 25.0 1.34
## 9 09 38.2 0.161 23.1 1.15
## 10 10 38.1 0.209 20.8 1.16
## 11 11 38.0 0.200 18.1 1.03
## 12 12 38.0 0.178 16.1 0.845Table of anomaly regression analyses (2000-2009) - 0m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano -0.00188507 0.003123799 0.54649521 3.779479 6.263081
## Temperature_ano 0.03810298 0.018994081 0.04541915 -76.394730 38.082303
## P int. F df R2 P value
## Salinity_ano 0.54649562 0.364157 474 0.0007676739 0.54649521
## Temperature_ano 0.04541937 4.024216 474 0.0084184354 0.04541915Plot of anomalies : Temperature (2000-2009) - 0m
Plot of anomalies : Salinity (2000-2009) - 0m
Jan.2010 - Dec.2022
Table of monthly means (2010-2022) - 0m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.0 0.375 14.5 0.711
## 2 02 38.0 0.233 13.7 0.375
## 3 03 37.9 0.245 13.7 0.482
## 4 04 37.8 0.295 15.1 0.927
## 5 05 37.8 0.254 17.7 1.59
## 6 06 37.7 0.296 21.8 2.02
## 7 07 37.8 0.373 25.0 1.76
## 8 08 38.0 0.450 25.6 1.45
## 9 09 38.2 0.180 23.7 1.53
## 10 10 38.2 0.179 20.9 1.10
## 11 11 37.9 0.518 18.3 1.00
## 12 12 38.0 0.458 16.3 0.949Table of anomaly regression analyses (2010-2022) - 0m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano 0.003813676 0.003557838 2.841797e-01 -7.689935 7.174074
## Temperature_ano 0.092233433 0.013161962 6.362107e-12 -185.980424 26.539958
## P int. F df R2 P value
## Salinity_ano 2.841805e-01 1.148987 618 0.001855753 2.841797e-01
## Temperature_ano 6.362608e-12 49.106106 618 0.073610638 6.362107e-12Plot of anomalies : Temperature (2010-2022) - 0m
Plot of anomalies : Salinity (2010-2022) - 0m
DEPTH 50
Jan.2000 - Dec.2009
Table of monthly means (2000-2009) - 50m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.1 0.0914 14.2 0.654
## 2 02 38.1 0.0779 13.5 0.461
## 3 03 38.1 0.0799 13.4 0.473
## 4 04 38.1 0.0846 13.8 0.528
## 5 05 38.0 0.118 14.8 0.859
## 6 06 38.0 0.111 14.9 0.622
## 7 07 38.0 0.113 14.9 0.613
## 8 08 38.0 0.137 15.2 0.528
## 9 09 38.0 0.102 15.8 1.10
## 10 10 38.1 0.176 18.3 1.92
## 11 11 38.1 0.144 18.0 1.02
## 12 12 38.0 0.131 16.0 0.813Table of anomaly regression analyses (2000-2009) - 50m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano -0.002918298 0.001847716 0.11490525 5.851096 3.704616
## Temperature_ano 0.027183103 0.013892758 0.05097301 -54.501262 27.854567
## P int. F df R2 P value
## Salinity_ano 0.11490562 2.494531 477 0.005202419 0.11490525
## Temperature_ano 0.05097325 3.828434 477 0.007962162 0.05097301Plot of anomalies : Temperature (2000-2009) - 50m
Plot of anomalies : Salinity (2000-2009) - 50m
Jan.2010 - Dec.2022
Table of monthly means (2010-2022) - 50m data
## # A tibble: 12 × 5
## month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.1 0.134 14.5 0.663
## 2 02 38.1 0.128 13.8 0.356
## 3 03 38.1 0.146 13.6 0.349
## 4 04 38.0 0.124 14.2 0.517
## 5 05 38.1 0.116 15.0 0.703
## 6 06 38.1 0.0869 15.1 0.607
## 7 07 38.1 0.0752 15.3 0.643
## 8 08 38.0 0.0593 15.4 0.640
## 9 09 38.1 0.0984 16.2 1.92
## 10 10 38.1 0.152 18.9 1.81
## 11 11 NA NA 18.0 1.19
## 12 12 NA NA 16.3 1.00Table of anomaly regression analyses (2010-2022) - 50m data
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano 0.0114991 0.001276746 3.963691e-18 -23.18538 2.574279
## Temperature_ano 0.0257115 0.010697870 1.653841e-02 -51.84285 21.570465
## P int. F df R2 P value
## Salinity_ano 3.964160e-18 81.118272 523 0.134275482 3.963691e-18
## Temperature_ano 1.653859e-02 5.776438 613 0.009335259 1.653841e-02Plot of anomalies : Temperature (2010-2022) - 50m
Plot of anomalies : Salinity (2010-2022) - 50m
Comparison SOMLIT - EOL results : 2018-2022
SURFACE
temperature (°C) & salinity
Table and plot of monthly means (2017-2022) - EOL temperature & salinity
monthly_means_EOL_1822 <- ungroup(EOL_hourly_1822) %>%
group_by(Month) %>%
summarise(
Salinity_month = mean(salinity, na.rm = TRUE),
SD_sal = sd(salinity, na.rm = TRUE),
Temperature_month = mean(temperature, na.rm = TRUE),
SD_temp = sd(temperature, na.rm = TRUE))
monthly_means_EOL_1822## # A tibble: 12 × 5
## Month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.0 0.107 14.6 0.564
## 2 02 38.1 0.151 13.9 0.277
## 3 03 38.0 0.268 13.9 0.388
## 4 04 37.9 0.303 15.1 1.08
## 5 05 37.8 0.202 18.2 1.92
## 6 06 37.8 0.231 22.4 1.80
## 7 07 37.9 0.150 26.1 1.40
## 8 08 38.1 0.195 26.7 1.01
## 9 09 38.2 0.167 24.5 1.44
## 10 10 38.1 0.182 21.2 1.10
## 11 11 37.7 3.13 18.7 1.05
## 12 12 38.0 0.160 16.4 0.765monthly_means_somlit_1822 <- ungroup(SOMLIT_1822) %>%
group_by(Month) %>%
summarise(
Salinity_month = mean(impute_sal, na.rm = TRUE),
SD_sal = sd(impute_sal, na.rm = TRUE),
Temperature_month = mean(impute_temp, na.rm = TRUE),
SD_temp = sd(impute_temp, na.rm = TRUE))
monthly_means_somlit_1822## # A tibble: 12 × 5
## Month Salinity_month SD_sal Temperature_month SD_temp
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 01 38.1 0.118 14.6 0.600
## 2 02 38.1 0.124 13.8 0.310
## 3 03 38.0 0.194 13.8 0.405
## 4 04 38.0 0.151 14.8 0.718
## 5 05 37.7 0.296 17.9 1.91
## 6 06 37.5 0.548 22.4 1.88
## 7 07 37.6 0.426 26.1 1.25
## 8 08 37.8 0.556 26.2 1.51
## 9 09 38.1 0.220 24.4 1.35
## 10 10 38.1 0.238 21.2 1.06
## 11 11 38.0 0.418 18.9 1.02
## 12 12 38.0 0.191 16.3 0.754Table of anomalies regression analyses (2018-2022) - EOL temperature & salinity
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano 0.002164672 0.003343598 5.173718e-01 -4.373596 6.755563
## Temperature_ano 0.031911638 0.004104452 7.733525e-15 -64.477267 8.293022
## P int. F df R2 P value
## Salinity_ano 5.173726e-01 0.419137 38850 1.078848e-05 5.173718e-01
## Temperature_ano 7.733608e-15 60.448817 39611 1.523736e-03 7.733525e-15Table of anomalies regression analyses (2018-2022) - SOMLIT temperature & salinity
## Slope SE Slope P Slope Intercept SE int.
## Salinity_ano -0.004897158 0.01386080 0.7241470 9.894662 28.00563
## Temperature_ano 0.039645714 0.04992486 0.4278683 -80.103804 100.87275
## P int. F df R2 P value
## Salinity_ano 0.7241470 0.1248278 256 0.000487371 0.7241470
## Temperature_ano 0.4278684 0.6306070 256 0.002457255 0.4278683EOL : Plot of anomalies : Temperature (°C) (2018-2022)
EOL : Plot of anomalies : Salinity (2018-2022)
without interpolation
SOMLIT : Plot of anomalies : Temperature (°C) (2018-2022)
SOMLIT : Plot of anomalies : Salinity (2018-2022)
